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ABSTRACT 


The variable-interval time-averaging (VITA) technique developed by 
Blackwelder and Kaplan [J. Fluid Mech. 7b, 89 (1976)] is applied to 
data obtained from large-eddy simulation of turbulent channel flow in 
an investigation of the organized structures associated with the burst- 
ing phenomenon in the near-wall region. Conditionally averaged veloci- 
ties, shear stress, pressure, and vorticity are discussed in conjunc- 
tion with the bursting phenomenon detected by the VITA technique. The 
conditionally averaged pressure reveals that the ejection process is 
associated with a localized adverse pressure gradient. In the plane 
perpendicular to the flow direction, the conditionally averaged vortic- 
ity field indicates that a pair of counterrotating streamwise vorticity 
is being lifted through the ejection process. 
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I. INTRODUCTION 


During the past decade, coherent structures in turbulent flow have 
been the subject of considerable interest among turbulence researchers. 
For a complete review of these structures, the reader should refer to 
the recent review article by Cantwell.^ In the present paper, the 
coherent structures associated with wall-bounded turbulent flows in the 
near-wall region are investigated. Substantial experimental evidence 
has been collected which indicates that there exist quasi-cyclic, 
spatially coherent structures in the vicinity of the wall in turbulent 
boundary layers. Although there is little disagreement about the 
existence of such organized structures, the lack of consensus on the 
detailed description of the structure still prevails (for example, see 
Kline and Falco^). 

The bursting phenomenon is a subject of intense interest to 
researchers who are working to identify the organized structures in the 
wall-bounded shear flows. It is composed of a sequence of quasi-cyclic 
events that occur in the wall region. Extensive experimental work on 
the bursting phenomenon has been reported in the literature: e.g. , 

Kline et al.,^ Corino and Brodkey,** Brodkey et al.,® and Of fen and 
Kline.® Blackwelder and Kaplan^ used a novel technique, called VITA 
(variable-interval time-averaging) , to single out the bursting event 
from the usual chaotic turbulent flow. They then obtained the condi- 
tionally averaged statistics and reported some organized structures 
associated with the bursting phenomenon. In the present study, the 
same technique was applied to the data obtained from the large-eddy 
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simulation of turbulent channel flow (Moin and Kim®), and the results 
are compared with the experimental results. In addition, the informa- 
tion on pressure and vorticity; which is not available from laboratory 
experiments, is analyzed in detail to gain a better understanding of 
the bursting phenomenon. In the aforementioned work of Moin and Kim, 
a fully developed turbulent channel flow was simulated numerically at a 
Reynolds number of 13,800, based on the centerline velocity and channel 
half-width. The computed flow field was used to study the statistical 
properties of the flow as well as its time-dependent features. The 
agreement between the computed statistics and detailed flow structures 
and the experimental data was good. In addition, the turbulence struc- 
tures in the vicinity of the wall of the channel flow were found to be 
the same as the structures in the wall region of turbulent boundary 
layers. 

The terms "organized" and "coherent" structures have been used 
widely to imply different things by different investigators. It might 
be appropriate here to clarify their meaning in the present context to 
avoid a possible confusion. In the present paper, the words organized 
and coherent are used to imply structures that are recognizable through 
the ensemble-averaging process defined by the VITA technique. In other 
words, it is the structure that is ordered enough to be detected by 
the VITA technique. 

In Sec. II the VITA technique is described briefly, and the results 
of the conditionally averaged field are presented in Sec. III. A short 
discussion and a conclusion based on the present Investigation are 
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presented in Sec. IV and a short summary of the present study is pro- 
vided in Sec. V. 

II. CONDITIONAL SAMPLING 

The detection of the bursting phenomenon can be very subjective 
depending on the investigator's detection criterion. Most of the 
information about the phenomenon is descriptive and is primarily based 
on flow-visualization experiments. If one tries to identify the event 
quantitatively, subjective criteria, such as the choice of detection 
function and threshold values, have to be introduced. Different 
investigators have used different techniques and have produced different 
results. In the present study, the VITA technique was used because it 
was simple and easy to Implement, although not necessarily the best 
method available. The technique has been used previously by several 
investigators (Blackwelder and Kaplan,^ Blackwelder and Eckelmann,® and 
Antonia et al. to name a few). The results of the present study can 
be compared with those of the previous ones. A brief description of 
the technique is given below, but the reader should refer to Ref. 7 for 
details. 

The variable-interval time-average of a fluctuating quantity 
Q(x^,t) is defined as 

J t+T/2 

Q(x^,T)dT , (1) 

c-T/z 

where T is the averaging time. Note that 

lim Q = Q , (2) 

T-xo 
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where the bar indicates the conventional time-average. If one wants 
to obtain a local average of some particular phenomenon, the averaging 
time T must be of the order of the time-scale of the phenomenon under 
study. To represent a measure of the turbulent energy during time T, 
a localized variance is defined in the following way: 

■^(x^,t,T) = u^(x^,t,T) - [u(x^,t,T)]^ , (3) 

where u is the fluctuating component of the streamwise velocity. The 
detection criterion for the bursting phenomenon is then completed by 
using a threshold level on the VITA variance. The detection function 

~v 

D(t) is defined as ^ 

( 1, for v^ > k uf_ and 9u/9t > 0 , 

( 4 ) 

0, otherwise 

where k is the threshold level and u is the root mean square of 

rms 

the fluctuating streamwise velocity. The second constraint, 9u/9t > 0, 
was not a part of the original detection function used by Blackwelder 

•j 

and Kaplan, but was added in the subsequent work of Chen and 
Balckwelder. 

The conditional average of a quantity Q is defined by 

N 

(Q(x^.T))y+ =1^ QCx^.tj + T) , (5) 

where x^ is the position at which the sampling occurred, the subscript 
y*^ indicates the position at which detection occurred, and N is the 
total number of the events detected by Eq. (4). 

To obtain a spatial structure rather than the temporal structure, 
in the present study, the conditional averaging process was modified 
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in the following way to give the variable-interval space-averaging: 


x+L/2 


Q(x 


I u/ 

5y>ZjtQ,L) — ^ I Q (C »y » z , tQ)d^ , 


'x-L/z 


2 ^2 
var = u - u 


D(x) = 


=r- 

10. 


for var > k u and 8u/8x < 0 

ms 

otherwise 


and 


( 6 ) 

(7) 

(8) 


N 

(Q(?.y.z,to))y+ = ^ ^ + 5,y,z,to) , (9) 

where L is the width of the spatial averaging. (Note that most of 
the quantitative measurements produce temporal structures and that 
flow-visualization experiments generally yield spatial structures.) 

For all the data discussed in this paper, a threshold value of 1.2 and 
L = 8Ax, corresponding to 500 wall units (i.e., L"^ = Lu^/v = 500), were 
used. The location of the detection probe was set at y"^ = 21. These 
values are comparable with the values used by the previous investigators. 
Blackwelder and Kaplan,^ for example, used k = 1.2 and T"*" E Tu^/v = 10 
with the detection probe at y"^ = 15. In addition, the ensemble- 
averaging process was taken over several flow fields at different time 
to obtain better statistics. Most of the results reported here were 
results of the ensemble averages over more than 1000 independent events. 


III. CONDITIONAL AVERAGES 


Figure 1 shows the conditionally averaged velocity profiles as 
functions of y"^ and the streamwise locations relative to the point of 
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the detection. The location of the detection corresponds to 5=0. 

The negative 5 indicates upstream points of the detection, and the 
positive 5 indicates downstream points of the detection. Each incre- 
ment of 5 corresponds to about 62 wall units (i.e., = 62).; The 

dashed lines in the figure represent the mean-velocity profile and the 
symbols are the conditionally averaged profiles around the event. 

Recall that in contrast to the results of Blackwelder and Kaplan^ (BK, 
in short, hereinafter), the current results display the spatial struc- 
tures of the event rather than the temporal structures. One can relate 
the present results to those of BK using the Taylor’s hypothesis, in 
which case, the downstream separations correspond to the negative time 
delays and the upstream separation to the positive time delays; that is, 
one should scan Fig. 1 from the bottom right to the top left. It 
appears the agreement between the two results is good. Both results 
show clearly the sweep-ejection process (ejection-sweep in the case of 
BK) through the detection point. These conditionally averaged profiles 
are shown as functions of u"^ and y"^ in Fig. 2. In Fig. 2(a), the time- 
averaged wall-shear velocity u^, was used to obtain the nondimensional 
variables u"*" and y"^. With this nondimensionalization, the excess and 
the defect in the velocity during the event are confined within y"*" 
between 10 and 70. The outer parts of the profiles are not affected at 
all. In Fig. 2(b), the conditionally averaged wall-shear velocity, 

<u^>, is used to nondimensionalize u and y. With this nondimensional- 
ization, the conditionally averaged velocity profiles follow closely 
the mean-velocity profile up to about y"^ = 15, suggesting that the 
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velocity near the wall has adjusted rather quickly with the change in 
the wall shear. The deviation now is more pronounced in the outer part 
of the profiles. 

Conditionally averaged profiles of u, v, w, and uv are shown in 
Fig. 3. The total length in the streamwise direction in the figure is 
about 4000 wall units. Other than the fact that one has to follow the 
figures from the right to the left, these figures are similar to those 
of BK. The only difference between the two results is that the magni- 
tudes of the present results are consistently slightly lower than the 
experimental results. This could be attributed to the differences in 
parameters employed in the sampling process or, simply, to the differ- 
ence between the temporal and spatial coherence. In any event, the 
difference does not seem to be significant. Note that in Fig. 3(a), 
the bottom curve represents the streamwise velocity at y"*" = 1.8. 
Although it is not shown here, the conditionally averaged wall shear is 
very similar to this curve. One interesting feature of the curve is 
that it has a long tail to the upstream side (sweep side) . This causes 
the deviation from the mean to become more noticeable further along 
the upstream side in Fig. 2(b). 

The conditionally averaged pressure is shown in Fig. 4 together 
with the streamwise velocity. The pressure signatures show localized 
peaks associated with the event. As a result of these peaks, localized 
adverse pressure gradients are formed followed downstream by favorable 
pressure gradients around the detection point. Note that these peaks 
are more or less aligned in the direction normal to the wall. In con- 
trast, the peaks of u are skewed such that they arrive first away from 
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the wall. The skewness of the u-signature can be interpreted as 
follows: (1) either the large-scale structure moves toward the wall at 

an oblique angle; or (2) the convection velocity of the large-scale 
structures is a function of the distance from the wall. The pressure 
profiles, however, indicate that the convection velocity of the pres- 
sure signature associated with the large-scale structure is constant 
across the layer. This produces a phase lag in the velocity signature 
relative to the pressure signature near the wall. From Fig. 4, it 
appears that the ejection process (the slowing-down of the streamwise 
velocity in the figure) is associated with the localized adverse pres- 
sure gradient near the wall. Based on this observation, the following 
statements can be made to describe the bursting event detected by the 
VITA technique: (1) the sweeping motion of the large-scale structure 

moves toward the wall at an oblique angle; (2) the pressure wave asso- 
ciated with the sweep has localized peaks and precedes the velocity 
near the wall; (3) the fluid near the wall encounters the adverse pres- 
sure gradient and is ejected upward; and (4) as a result of this lift-up 
of the low-momentum fluid, the velocities away from the wall are 
decreased. Note that this description is consistent with the flow model 
of Of fen and Kline,® which describes the lift-up process as a response 
to a temporary local adverse pressure gradient. 

The conditionally averaged spanwise vorticity is shown in Fig. 5. 

In the regions shown in the figure, most contributions to the vorticity 
derive from -9u/9y, and the contribution of 9v/9x is negligible. 

The profiles show the excess and defect vorticity associated with the 
sweep-ejection motions. The corresponding streamwise and normal 
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vorticity were also computed; they showed no coherent motions in this 
C - y plane, as expected by symmetry. 

So far, we have examined the conditional averaged quantities in 
5 - y plane only. One can also study these profiles in 5 - z and 
y - z planes. In the 5 - z plane, we can study the spanwise extent 
of the structure associated with the event. In Fig. 6, the condition- 
ally averaged streamwise velocities, shear stresses, and pressure in 
this plane are shown at y"*" = 21 (note that this is the y-location of 
the detection point). The spanwise correlations of u and uv vanish 
beyond z"*” = ±30, whereas those of the pressure seem to extend to about 
±60. (The scaling of the large-scale structure was not captured prop- 
erly in the simulation of Moin and Kim® because of the numerical reso- 
lution problem; therefore, any quantitative information must be taken 
with some care.) This larger spanwise extent of the pressure eddies 
was also apparent from the contour plots of the instantaneous pressure 
(see Moin and Kim®). The conditionally averaged spanwise velocity 
profiles are shown in Fig. 7. The direction is such that the positive 
w represents moving downward (l.e., to the positive z) , and the nega- 
tive corresponds to moving upward (i.e., to the negative z). The 
profiles indicate the spanwise outward motions during the sweep and the 
inward motions during the ejection. These spanwise outward motions 
during the sweep, referred to as "splatting" by Moin and Kim,® resulted 
in a peculiar intercomponent energy transfer in the wall region. 

The conditionally averaged streamwise vorticity in y - z planes 
is shown in Fig. 8 at various streamwise locations. Note that flow 
goes into the plane in these figures; we then can investigate the 
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streamwise vorticity which is associated with the sweep-ejection pro- 
cess. The total spanwise width is about 2000 wall units. The sense of 
rotation is such that positive vorticity indicates clockwise rotation 
and negative vorticity indicates counterclockwise rotation. Figure 8(a) 
shows the profiles of the vorticity at C = -10, which is about 620 wall 
units upstream of the detection point (recall that C = 0 is at the 
center of the event and each increment corresponds to 62 wall units). 

At this upstream location, no detectable correlation is found. Fig- 
ure 8(b) is the vorticity signatures at ^ = -7, and we begin to see 
that an appreciable amount of the streamwise vorticity components sur- 
vived the ensemble-averaging process. 

Figures 8(c) through 8(k) show how the streamwise vorticity signa- 
tures change their patterns as we move downstream through the center of 
the event. Several points need to be addressed explicitly here. First, 
the streamwise vorticity appears as a counterrotating pair, and there 
seems to be two distinct groups: one near the wall and the other 

slightly away from the wall. These two groups have opposite signs of 
vorticity as a pair. Whether this streamwise vorticity is associated 
with the actual revolving vortex motion cannot be determined from the 
current figures. However, in the previous work of Moin and Kim,® it 
was found that the streamwise vorticity very close to the wall, say 
y"^ less than 10, was due to the gradient of spanwise velocity, whereas 
some of the vorticity at y"^ between 20 to 100 was Indeed associated 
with the revolving vortex motion. Next, the magnitudes of the vorticity 
near the wall Increase as we approach the detection point, and the 
vorticity that existed away from the wall slowly disappears; in fact at 
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C = 0, only the former group is noticeable. Also we note that the 
center of the pair of the vorticity (the y location of the maximum 
strength) of the former group starts at the wall at C = -5 and moves 
upward through the detection point (y*^ =21 is the center at this 
point) and moves further upward (y"^ =46 at ^ = 2) before it loses 
its identity far downstream. These figures suggest that the pair of 
counterrotating streamwise vorticity which is generated near the wall 
is being lifted-up through the ejection process. It is possible that 
this pair of counterrotating vorticity is in fact the two legs of the 
horseshoe vortex referred to by Offen and Kline^ in their flow model of 
the near-wall region, or the hairpin- type vortices observed farther out 
in the boundary layer by Head and Bandyopadhyay. The angle of the 
lifting vorticity with respect to the wall is roughly estimated to be 
about 10° at the center of the event, where the center of the vorticity 
is at y"^ = 21. This angle is far from the 45° of the hairpin vortices 
measured by Head and Bandyopadhyay . ^ ^ However, recall that we are 
still very close to the wall and that this angle is in agreement with 
the angle of large-eddy structures close to the wall reported by 
Rajagopalan and Antonia. 

Finally, one more point needs to be mentioned in conjunction with 
the figures. Downstream of the center of the event, a new pair of 
counterrotating vorticity appears, with the opposite sign of rotation 
relative to the pair of vorticity being lifted. This vorticity is the 
result of the upward motion and the no-slip boundary condition. 
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IV. DISCUSSION 


Lately some investigators have suggested that a turbulent boundary 
layer is composed almost exclusively of hairpin or horseshoe vortices. 
However, this may be too simplified a picture of the complicated turbu- 
lent flow. There exists some evidence that the hairpin and horshoe 
vortices do exist in the turbulent boundary layer, but they may not be 
the only structures that exist in the flow. Moin^** recently examined 
carefully the data of Moin and Kim® and also found some evidence of the 
hairpin-type vortices. However, he also found many other structures in 
addition to the motion of the hairpin-type vortices. The present study 
seems to suggest that the pair of counterrotating vorticity is a domi- 
nant structure during the bursting event (in the sense that they survive 
the ensemble-averaging process). Nevertheless, the possibility of 
whether this streamwise vorticity is indeed a part of the two legs of 
the hairpin- type vortices is still an open question. 

Based on the results presented in the previous section, one can 
make the following conclusions with regard to the bursting event 
detected by the VITA technique (a simple sketch characterizing this flow 
pattern is shown in Fig. 9): 

s 

1. High-speed fluid moves toward the wall at an oblique angle 
during the sweep. Unfortunately the cause of this sweeping motion is 
not known at present. 

2. The impingement of the high-speed fluid at the wall (splatting) 
and the viscous boundary condition create a pair of counterrotating, 
streamwise vorticity. 
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3. As this vortical fluid moves downstream, it encounters a 
localized adverse pressure gradient, which was set up by the sweeping 
motion (the pressure wave moves ahead of the sweeping fluid), and is 
ejected upward. 

4. As this pair of the streamwise vorticity is being lifted, they 
become a pair of counterrotating streamwise vortices (this notion of 
vortices rather than the vorticity, however, is inferred from the pre- 
vious study of Moin and Kim®). 

5. This pair of the counterrotating vortices scoops up more low- 
speed fluid in between the vortices, and the flow in the large portion 
of the inner layer decelerates. They also create a pair of streamwise 
vorticity of opposite signs with respect to them near the wall because 
of the viscous boundary condition. 

Note that this description of the bursting event is generally con- 
sistent with the flow model proposed by Blackwelder and Eckelmann® on 
the role of the streamwise vortices for the event. However, the stream- 
wise extent of these vortices is shorter than the one their model sug- 
gests. It appears that these streamwise vortices are much shorter than 
the wall layer streaks, suggesting that the streamwise vortices are 
needed only for initiating the formation of the streaks and that they 
do not have to be present side by side with the streaks. In addition, 
these streamwise vortices are tilted with respect to the wall. This is 
consistent with the result of Kim and Moin,^® which showed that the 
contour plots of the streamwise vorticity did not exhibit any streaky 
structure in contrast to the contours of the streamwise velocity. 
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V. SUMMARY 


This paper reports an attempt to identify some organized structures 
of wall-bounded shear flows using the results obtained from a three- 
dimensional time-dependent numerical simulation. The availability of 
the three-dimensional turbulent flow field allowed the author to inves- 
tigate the spatial structures associated with the bursting phenomenon. 
The VITA technique developed by Blackwelder and Kaplan^ was applied to 
the flow field to detect the bursting events. The present results 
agree well with available experimental data, if one converts the tem- 
poral structures of the experimental data to the spatial structure, 
using the Taylor's hypothesis. Upon this validation, the present study 
was extended to include conditionally averaged pressure and vorticity. 
The conditionally averaged pressure revealed that the ejection process 
at the end of the sweeping motion was associated with the localized 
adverse pressure gradient. This result was consistent with the flow 
model of Of fen and Kline,® although their model was based. on a specula- 
tion because of the lack of pressure data associated with their visual 
data. The current results confirmed their hypothesis that the lift-up 
of the inner layer may be akin to convected separation brought about 
by a temporary local adverse pressure gradient. 

The conditionally averaged profiles of streamwise vorticity in the 
y - z plane indicated that a pair of counterrotating vorticity was 
being lifted up through the ejection process. This streamwise vorticity 
was initially formed at the wall as the result of "splatting" and the 
viscous boundary condition. Although the present study indicated that 
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the pair of counterrotating streamwise vorticity is a dominant struc- 
ture for the bursting event, further work is required to determine 
whether this structure is related to the horseshoe vortex of Offen and 
Kline® or to the hairpin vortices of Head and Bandyopadhyay . 
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FIGURE CAPTIONS 


FIG. 1. Conditionally averaged and mean-velocity profiles with upstream 
and downstream separation relative to the point of detection. 

(ooo = <u>/U„; — , u/U„). 

FIG. 2. Conditionally averaged and mean-velocity profiles as functions 
of y"^ and u"^: (a) y"*" and u"*" are nondimensionalized by the mean-wall- 
shear velocity; (b) y"^ and u"*" are nondimensionalized by the condition- 
ally averaged wall-shear velocity (ooo = <u>/u^; — , = u/u^). 

FIG. 3. The conditional averages in C - y plane. The detection was 
applied at y"*" = 21: (a) streamwise velocities; (b) normal velocities; 

(c) spanwise velocities; (d) Reynolds shear stress. 

FIG. 4. The conditional averages of the streamwise velocity and pres- 
sure (- = <u>; — = <p>). 

FIG. 5. The conditional averages of the spanwise vorticity. 

FIG. 6. The conditional averages as a function of the spanwise coor- 
dinate at y"*" = 21: (a) streamwise velocity; (b) Reynolds shear stress; 

(c) pressure. 

FIG. 7. Conditionally averaged spanwise velocities as a function of 
the spanwise coordinate at y"*" = 21. 

FIG. 8. Conditionally averaged streamwise vorticity as functions of 
the normal and spanwise coordinates: (a) C = -10; (b) ^ = -7; (c) ^ = -5 

(d) 5 = -3; (e) K = -2; (f) 5 = -1; (g) C = 0; (h) ? = 1; (i) K = 2; 

(j) 5=3; and (k) 5=5. 
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FIG. 9. Flow patterns and associated vorticity during the sweep and 
ejection process: (a) end view; (b) side view. 
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